Fluctuation Effect in Superconductivity in an Electron Model with d-Wave Attaction 
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On the basis of i-matrix approximation, we study the superconductivity in the tight-binding model 
with d-wave attraction. The low-lying collective modes are considered as the predominant long- 
range fluctuations in the system. The Green's function is obtained as an analytic solution to a cubic 
equation. The superconducting order parameter and the transition temperature are substantially 
reduced from the values of the mean-field theory. The calculated phase boundary of the supercon- 
ductivity can reasonably describe the experiment results for the cuprates. The results for density of 
states and the pseudogap parameter are presented. 
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The fluctuation effect takes important role for describ- 
ing the superconductivity in a quasi-two-dimensional su- 
perconductor with low carrier density [1]. The super- 
conducting order parameter and thereby the transition 
temperature Tc can be substantially suppressed from the 
values given by the mean- field theory (MFT). In most 
cases, the long-range pairing fluctuation is important as 
it prohibits the ofF-diagonal-long-range ordering in sys- 
tems of dimensions < 2 [2] . For the under doped high- Tc 
cuprates (HTC), Emery and Kivelson attributed the sup- 
pression on Tc to the long-range classical phase fluctua- 
tion [3]. Above Tc, hole pairing becomes locally without 
long-range phase coherence. On the other hand, since 
there are preformed pairs above Tc, the superconductiv- 
ity below Tc can be viewed as a consequence of Bose- 
Einstein condensation [4-6] . Along with the approach by 
Leggett [4], much effort has been devoted to investiga- 
tion of the crossover from the weak-coupling BCS super- 
conductivity to the Bose-Einstein condensation (BEC) of 
bound pairs [7-16]. Most of the work has been performed 
for the s-wave pairing because of its computational sim- 
plicity. 

In this paper, we investigate the superconductivity in 
the tight-binding model with d-wave attraction. The 
low-lying collective modes of the boson fields above the 
condensate are treated as the predominant fluctuations. 
We present a self-consistent formalism for the Green's 
function based on the i-matrix approximation. The sin- 
gle particles and the collective modes are treated on the 
equal-footing manner. The new result is that below Tc, 
with the condensation taking place, the single pairs be- 
come moving collectively. Even at the ground state, there 
remains the zero-point motion. We calculate the phase 
boundary of superconductivity and compare the result 
with experiment data of HTC. We also study the pseu- 
dogap phenomena at Tc. 

The Hamiltonian of the electron system is given by 



where c|.^ {cka) is the creation (annihilation) opera- 
tor for electrons with momentum-fc and spin-a, ^k = 
~2t{coskx + cosfcy) — 2tz cos kz — /i with fi the chemi- 
cal potential, N the total number of lattice sites, Vkk' = 
—vrikTik' with rjk = cosk^ — cosky, and p{k,q) — 
C-k+q/2iCk+q/2] IS thc pair operator. The differences be- 
tween Eq. (1) and the t ~ J model are, (1) here only 
the d-wave attraction channel in the interaction is taken 
into account, (2) the prohibition of doubly occupation on 
the same site is released by assuming the hopping inte- 
grals t and tz are proportional to the hole concentration 
6, e.g., t — tgS with tg a constant. For the quasi- two- 
dimensional system, tz/t ^ 1 is supposed. This model 
has been adopted by a number investigators for studying 
the d-wave superconductivity as well as the pseudogap 
phenomena in cuprates [12-14,17]. 

In Nambu's representation, thc Green's function 
G(k,Zn) of thc electrons is given by 



G(fc, Z„) = [z„ - ^k<^3 - S(fc, Zn)Y 



(2) 



where z„ — iTTT{2n -f- 1), n is an integer, T the tempera- 
ture, and a the Pauli matrix. Throughout this paper, we 
use the units in which ?i = fc^ = 1- To express the self- 
energy, firstly, we note that the off-diagonal part comes 
from averaged boson fields of momentum q = 0. In the 
superconducting state, {p{k, 0)) is a macroscopic quantity 
as compared with any other pair fields else. Therefore, 
the predominant contributions are the static mean filed 



^Mk, Zn)^^Y. ""kk' {P{k', 0)) = Ak 



(3) 



koc kk' q 
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and its complex conjugate I]2i(fc,z„). For our uniform 
system, we suppose A^ is real. The quantity A^ = Aryfc 
should be differentiated from that of thc MFT since the 
fluctuation effect is under consideration in the present 
Green's function. Secondly, we take into account of the 
pair fluctuation terms g 7^ of the interaction in the 
diagonal part. By the t-matrix approximation, which is 



the infinite sum of ladder diagrams, the diagonal part of 
the self-energy is given by 



qrn 



where Zm = i2TrTm, m is an integer, /x = 1, 2 with 1 — 2 
and 2 = 1, and the t matrix W{q, Zm) is given by 

W{q, Z„,) ^v''[l + vxiq, ZmT'xiq, Z^) (5) 

where W and x s-re 2x2 matrices, with the elements of 
X defined by 



T 



XiJ-y 



u{q, Z,n) = -r^ X! VkGfi,,{k+q/2, z„+Z„)Gp^(A:-(7/2, z„) 



kn 



(6) 



The chemical potential ji is determined by 

2T 

— Y^G^^{k,z^)e^-^^l-5 (7) 

kn 

where rj is an infinitesimal positive number. These equa- 
tions (2)- (7) form the closed system that self-consistently 
determines the Green's function and the t matrix. 

Two points about the present formalism need to be 
emphasized. Firstly, the existence of Goldstone mode in 
the superconducting state requires that the pair suscep- 
tibility x((7, Zm) should satisfy the condition [18] 



det|l + i;x(0,0)| = 0. 



(8) 



This equation is exactly consistent with Eq. (3). Any 
improper treatment of off-diagonal self-energy leads to 
violation of this consistency. At Tc, equation (8) reduces 
to the Thouless criterion [19]. Secondly, because of the 
divergence of W^^((7, Zm) s.t q — Q and Zm = 0, the di- 
agonal self-energy takes into account of the predominant 
fluctuation effect. 

It is a tremendous task to numerically solve Eqs. (2)- 
(7) because many multi-dimensional integrals over mo- 
mentum and the summation over Matsubara's frequency 
need to be computed in each iteration. However, since 
W{q, Zm) is strongly peaked with a divergence about 
g = and Zm = 0, the diagonal self-energy can be ap- 
proximately given by [14,20] 

S^^(fc, zn) ~ -rilGf,-^{k, z„)^ Y. '^A<A<(9' Zm)e'^-^-^ 

qm 

(9) 
where ^ ' means the summation over q runs a small q 
region, and the convergent factor g"*'^™'' with ai — \ and 
a2 = — 1 has been introduced. This convergent factor 
comes from the fact that the Green's function Gpp(fc — 
q, Zn — Zm) in the summation in Eq. (4) is connected with 
the effective interaction W^^{q,Zm)- The summations 
over q and m in Eq. (9) give rise to a constant 



T 



r' = -^E'^MM(9,^m)' 



.a^Z^ri 



(10) 



At small q and Zm, the t matrix can be approximated by 
the collective modes. At T < T^ the collective modes are 
sound- like waves with energy fig oc g, while at T^ fig oc 
<7^, the excitations are single pairs. The constant F^ is 
essentially a measure of the density of these uncondensed 
pairs. F is named as pseudogap parameter since at T^ 
there remains a gap in the density of states (DOS) at the 
Fermi energy. The g-integral in Eq. (10) is over a cylinder 
region in the momentum space. Since Q.q depends weakly 
on the out-plane wave number q^ , the integral over qz can 
be taken in the range (— 7r,7r). The cutoff qc for the in- 
plane wave number is determined such that the largest 
in-plane energy Vtq^ = 2\fE? + T^ since the collective 
mode is meaningful only within the gap. 

Now, note that two equations from the diagonal parts 
of Eq.(2) with S^^ given by Eq. (9) form the closed sys- 
tem for determining the diagonal Green's functions. By 
eliminating one of them, we can obtain a cubic equation 
for Go or G3. By introducing a function j/(fc,z„), the 
Green's function is obtained as 

G{k,z^) = [z„ + 3Afeai/(l + y)+aT3](2-y)/3F2. (11) 

where F^ = Fry^. The function y{k, Zn) is a real root to 
the cubic equation 



y^ -3Py-2Q = 



(12) 



where P = 1 + 3(F2 - AD/ie^ - 4) and Q = 1 + |(F2 
2Al)/{(_l — z^). The explicit form of j/(fc, z„) reads 




2VPcos(^/3), 



Q - Vd, D>0 
D <0 



(13) 



arccos((5/-\ 



qm 



and if — a.TCCOs{Q / vP^) ■ The 



where D = Q^ - P^, 

boundary condition y ^ 2 + 3F|/(^^ — z^) at |z„| -^ 00 
is useful for analytic contituation to real frequency z„ -^ 
Lu + irj. 

For describing the cuprates, we take v/2to — 0.1 and 
tz/t ~ 0.01 [13]. For LazCuO, v ~ 0.13 eV has been 
determined by experiments [21]. Therefore, our choice of 
v/2tQ corresponds to to — 0-65 eV, which is consistent 
with estimates from experiment data [13,22]. The small 
quantity tz/t describes the interlayer weak coupling and 
gives rise to a ^-freedom energy in Qq. This weak cou- 
pling prevents the summation over q in Eq. (10) from 
a logarithm divergence at the q = limit and ensures a 
finite transition temperature Tc- 

The result for Tc as a function of hole concentration S 
(sohd hue) is plotted in Fig. 1. The MET resuh (with 
the same scale as the present theory) and the experiment 
data [23] are also shown for comparison. The maximum 
transition temperature T^^Max ~ 108 K obtained by the 
present theory appears at a certain 6 between 0.125 and 
0.15. The improvement of the present theory to the MET 
is significantly. The phase boundary of superconductivity 



given by the present theory is close to the experiment 
results. In the optimally to over doped region, the present 
theory fits the experiment data very well. Contrary to 
the MFT in the under doped region, Tc of the present 
calculation increases reasonably with 5. Such a behavior 
of Tc reflects the feature of BEC [19]. 

There is still obvious discrepancy between the present 
theory and the experiment at small 5. This may date 
from the crude treatment of the short-range pair correla- 
tions. Local pairing without long-range phase coherence 
is not fully taken into account in the present model. Be- 
sides this, the short-range antiferromagnetic coupling is 
not correctly counted in. To describe the antifcrromeg- 
netism in cuprates at very small (5, one needs to restart 
with the t — J model. In addition, long-range Coulomb 
interaction may also take effect if there is no adequate 
screening. To what extent all these factors affect the su- 
perconductivity may be a subject for future study. 

The reduction of Tc from its MFT value stems from 
that a part of the pairs occupy collective modes. As men- 
tioned above, the density of these pairs is measured by 
r^. In Fig. 2, the result for the pseudogap parameter F 
at Tc is compared with the MFT Amf{Tc) as well as the 
pseudogap energy Eg determined by experiments [24]. 
The experimental observations indicate that Eg slightly 
depends on T above Tc [25]. We here plot the result for 
Eg to view the overall magnitudes only. The parame- 
ter r similarly as Amf(T'c) decreases monotonously with 
6. The larger ratio T / /^mf{Tc) at smaller 5 implies the 
larger occupation of the uncondensed pairs. This is con- 
sistent with the more significant reduction of Tc in the 
under doped region. 

In Fig. 3, we show the numerical results for A and 
F as functions of T at 5 = 0.1. The superconducting 
gap opens below Tc and reaches its maximum at T = 0. 
The parameter F decreases with T decreasing and re- 
mains a finite value at T = 0. At the ground state, 
the fluctuation effect comes mainly from the zero-point 
motion (quantum fluctuation) of the pairs. Generally 
speaking, in a quasi-two-dimensional system, the order 
parameter of the broken-symmetry state can be consider- 
ably suppressed by the quantum fluctuation. This can be 
confirmed also by the perturbation calculations [26-28]. 
The Chicago group [12-14] has considered the BCS-BEC 
crossover problem. By their approach, the uncondensed 
pairs below Tc are single pairs with energy Q,q oc q^, 
and completely condensation takes place at T = 0. By 
the present treatment, however, the BEC from the single 
pairs begins to occur at Tc^ while below Tc, with the con- 
densation taking place, the uncondensed pairs become 
moving collectively. At T = 0, the zero-point motion 
prevents pairs from fully BEC. 

The pairing fluctuation effect can be clearly re- 
flected in the DOS of single quasiparticles, p{E) — 
— 2i^j, ImGii(fc, E + iri)/TTN. A direct effect is the pair- 
ing fluctuation introduces lifetime to the single quasipar- 



ticles. Shown in Fig. 4 are the results for the DOS at 
5 = 0.1 and T/Tc = 0.2 and 1. p{E) depends on E lin- 
early at small E. This is in the character of d-wave gap. 
In contrast to the well-known MFT, the peaks in the 
DOS are broadened even at T < Tc due to the collective 
modes. The width scales with F. At Tc, p{E) shows still 
the existence of the d-wave gap with magnitude about 
2F. That is the pseudogap. To understand better of 
the psedogap, we consider the spectral function at Tc, 






6)/2^r2 



Aik,E) = ^(e,+^n-E')/iE^ 

[20] which is nonzero only for E'^ - 4F^ < ^l < E'^, 
with the noninteracting delta function peak becoming a 
square root singularity. For E"^ — 4F^ < 0, a part of spec- 
tral weight is shifted away. This results in the formation 
of a pseudogap in the DOS. 

In summary, we have investigated the superconductiv- 
ity in the tight-binding model with d-wave attraction. 
On the basis of i-matrix approximation, we have devel- 
oped the self-consistent formalism for the Green's func- 
tion. The analytic Green's function is given by Eq. (13). 
The low-lying collective modes are treated as the predom- 
inant long-range pairing fluctuation in the self-energy. 
The pairing fluctuation results in lifetime effect for the 
single particles below Tc and a pseudogap in the density 
of states at Tc- The transition temperature is substan- 
tially suppressed from its mean-field value. The phase 
boundary of superconductivity given by the present the- 
ory is close to the experiment results for the cuprates. 
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FIG. 1. Tc as a function of S. The solid and dashed 
lines represent the results of present approach and the 
MFT, respectively. The symbols indicate the experiment 
data for cuprates [23]: Yi_xCaxBa2Cu306 (solid squares), 
Yo.9Cao.iBa2Cu307-y (open squares), La2-xSrxCu04 (open 
diamonds), Yi_xCaxBa2Cu306.96 (solid diamonds), and 
YBa2Cu307-y (open circles). 
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FIG. 2. r (solid line) and Amf (dashed line, by the MFT) 
as functions of f) Pit TT,. The svmhols renresent the nsendop'an 
energT 




FIG. 3. A and P as functions of T at 5 = 0.1. Amf given 
by th 




FIG. 4. Density of states p{E) at T/Tc = 0.2 and 1. The 
hole concentration is S = 0.1. 



